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2-D  Design  of  Schottky  Diodes 

Jesus  Grajal,  Daniel  Moreno  and  Viktor  Krozer 


Abstract — The  optimization  of  Schottky  diode  operation 
together  with  the  design  of  circuits  based  on  these  devices 
can  be  accomplished  with  a developed  2-D  drift-diffusion 
simulator.  Evaluation  of  the  impact  of  different  geometries 
on  the  device  performance  and  study  of  limiting  processes 
like  velocity  saturation  at  high  frequencies  are  analysed  with 
this  simulation  tool. 


I.  Introduction 


Schottky  Contact  . 
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GaAs  Schottky  barrier  diodes  continue  being  the  key  ele- 
ments for  mixers  (heterodyne  receivers)  and  frequency  mul- 
tipliers (all-solid-state  power  sources)  for  millimeter  and 
sub-millimeter  bands.  The  key  points  in  the  progress  of 
the  performance  of  circuits  based  on  Schottky  diodes  have 
been  the  enhanced  physical  insight  into  and  the  optimi- 
zation of  submm-wave  Schottky  diode  operation  [1],  the 
improvement  in  frequency  multiplier  and  mixer  analysis 
methods,  and  in  physical  analytical  Schottky  diode  models 
[2],  [3],  as  well  as  numerical  physical  device  models  [4].  The 
last  milestone  in  this  progress  is  the  coupling  of  numerical 
physical  models  and  circuit  simulators  which  enables  the 
co-design  of  circuits  and  devices  [4],  [5].  These  advanced 
tools  can  only  be  effective  with  accurate  physical  models 
which  allow  the  analysis  of  different  structures  without  em- 
pirical adjustment  parameters.  Furthermore,  these  models 
must  account  for  the  limiting  mechanisms  at  high  frequency 
and  high  power  operation.  Besides,  these  simulators  must 
be  able  to  evaluate  the  impact  of  different  geometries  on 
the  device  performance.  This  is  not  possible  using  an  1-D 
simulator  developed  in  the  past  [5]. 

We  have  developed  a 2-D  numerical  simulator  based 
on  the  drift-diffusion  theory  which  incorporates  the  accu- 
rate physical  model  considered  in  the  1-D  numerical  simu- 
lator: impact-ionization,  non-constant  recombination  ve- 
locity, self-consistent  incorporation  of  the  tunnelling  and 
image-force  effects.  The  geometrical  aspects  are  included 
by  the  two-dimensional  nature  of  the  simulations.  Figure  1 
shows  the  different  structures  (half  of  the  device  exploiting 
axial  or  plane  symmetry)  which  can  be  analysed  with  this 
simulator.  It  is  also  clear  that  it  is  possible  to  study  pla- 
nar diodes  with  the  ohmic  contact  at  the  top  of  the  device 
(ohmic  contact  identified  with  symbol  (0),  figure  1),  and 
Whisker  and  quasi-vertical  planar  diodes  with  the  ohmic 
contact  at  the  rear  (ohmic  contact  identified  with  symbol 
(1),  figure  1). 

This  simulator  was  validated  by  comparison  of  simulated 
device  characteristics  with  experimental  results  obtained 
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Fig.  1.  Possible  geometries  in  the  2 D simulator,  (z)  identifies  the 
position  of  the  ohmic  contacts. 

for  submm-wave  Schottky  diodes  fabricated  at  the  Univer- 
sity of  Darmstadt  ( TUD ).  After  this  initial  step,  we  fo- 
cused on  the  impact  of  different  geometries  on  the  device 
electrical  characteristics  and  also  on  the  study  of  some  li- 
miting physical  processes  at  high  frequencies  like  velocity 
saturation.  The  objective  of  this  paper  is  to  show  that 
this  simulator  is  an  appropriate  tool  for  both  the  study  of 
physical  properties  of  these  devices  and  the  design  of  real 
devices. 

TABLE  I 

Parameters  for  Schottky  diodes  from  TUD  and  UVa. 


Diode 

Epitaxial  Layer 

Vbd)DC  (V) 

doping  (cm  6) 

length  (fim) 

D1038 

3 x 1017 

0.560 

-8 

UVa  6P4 

3.5  x 10ie 

1.0 

-20 

D1038  subs,  doping:  2 x 1018cm  3 length:  80-100  pm 


II.  Geometrical  aspect  of  the  design  of 
Schottky  diodes 

Impact  of  the  geometry  on  the  series  resistance  Rs  and 
the  capacitance  C of  Schottky  diodes  can  be  directly  stu- 
died with  this  simulator.  These  parameters  are  important 
because  their  optimization  can  improve  the  frequency  per- 
formance of  these  devices,  which  is  often  measured  by  the 
cut-off  frecuency:  fc0  = where  Cjo  = C(V&ias  = 

0). 

Figure  2 shows  the  series  resistance  of  a cylindrical 
Schottky  diode,  D1038 , whose  doping  profile  is  described 
in  table  I.  The  ohmic  contact  is  at  the  rear  of  the  de- 
vice and  the  radius  of  the  Schottky  contact  is  rsc  = 5 ixm 
(0  = 2-  rsc  in  figure  1)  for  different  substrate  lengths  and 
radii  of  the  ohmic  contact  ( r0hmic  = c + rsc  in.  figure  1).  It 
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is  important  to  notice  that  the  distribution  of  the  current 
in  the  substrate  changes  depending  on  the  relationship  bet- 
ween h (/i  & L substrate')  and  t ohmic  (figure  1),  and  this  is 
the  origin  in  variations  of  the  series  resistance.  For  short 
substrates  the  resistance  does  not  depend  on  the  radius 
of  the  ohmic  contact.  We  have  also  concluded  that  the 
resistance  does  not  depend  on  the  size  of  ohmic  contact, 
regardless  the  length  of  the  substrate  for  c > (3  - 4)  * rsc • 
Similar  results  apply  for  for  another  Schottky  diode  with 
rsc  = 2.5  fim  as  can  be  observed  in  figure  3. 
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Fig.  2.  Series  resistance  vs.  radius  for  the  ohmic  contact  for  different 
substrate  lenghts.  The  ohmic  contact  is  at  the  rear  of  the  device 
and  radius  of  the  Schottky  contact  is  rsc  = f = 5 /xra. 
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Fig.  3.  Series  resistance  vs.  radius  for  the  ohmic  contact  for  different 
substrate  lenghts.  The  ohmic  contact  is  at  the  rear  of  the  device 
and  radius  of  the  Schottky  contact  is  rsc  = f = 2.5  pm. 

A similar  study  has  been  performed  for  this  device  with 
the  ohmic  contact  at  the  top  for  different  substrate  leng- 
hts and  distances  between  the  contacts:  c — 3 jum,  as  can 
be  seen  figure  4.  For  this  structure,  the  series  resistance 
decreases  for  increasing  substrate  lengths  because  the  area 


for  the  current  flux  increases.  The  differences  among  the 
analysed  structures  are  nearly  negligible  for  distances  bet- 
ween contacts  less  than  5 jim. 
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Fig.  4.  Series  resistance  vs.  distance  between  contacts  (c  - 3 jwm, 
figure  1)  for  different  substrate  lenghts.  The  ohmic  contact  is  at 
the  top  of  the  device  with  radius  r0hmic  = 3 fim.  The  radius  of 
the  Schottky  contact  is  rsc  = f = 2.5  fim. 

The  capacitance  of  a small-area  submillimeter-wave 
diode  is  strongly  affected  by  the  edge  effects  as  has  been 
pointed  out  by  several  researchers  [6].  Figure  5 presents 
the  correction  factor  7 c for  the  capacitance  due  to  edge 
effects  [6]  as  a function  of  the  radius  of  the  anode  for  seve- 
ral doping  densities  in  the  epitaxial  layer.  It  is  clear  that 
the  edge  effect  only  impacts  on  the  characteristics  of  the 
device  if  the  radius  is  small  enough  compared  to  the  length 
of  the  epitaxial  layer  Wepi  (ffib  < 3 — 4),  which  is  typical 
for  submillimeter  varactors.  Of  course,  the  influence  of  the 
edge  effects  also  depends  on  the  doping  density  of  the  diode 
as  can  be  observed  in  figure  5. 

The  accurate  calculation  of  the  capacitance  is  fundamen- 
tal for  the  design  of  varactors  working  as  frequency  multi- 
pliers, where  the  basic  non-linearity  for  generating  harmo- 
nics is  the  capacitance.  Besides,  the  design  of  the  device 
geometry  in  order  to  reduce  the  series  resistance  will  im- 
pact dramatically  on  the  efficiency  of  a multiplier. 

III.  High  frequency  limiting  mechanisms: 
Velocity  saturation 

Besides  geometrical  aspects,  a numerical  simulator  is  also 
a useful  tool  to  gain  physical  insight  into  the  diode  opera- 
tion. High  frequency  operation  of  Schottky  diodes  working 
as  frequency  multipliers  or  mixers  is  limited  by  the  velo- 
city saturation  effect  [1].  It  has  been  outlined  in  [7]  that 
the  current  saturation  mechanism  should  not  occur  for  well 
designed  diodes  at  all  operating  frequencies.  This  can  be 
achieved  by  choosing  a sufficiently  high  doping  density  in 
the  device  No  > Jdisp/QVmax  with  Jdisp  being  the  displace- 
ment current  and  vmax  the  maximum  velocity  of  electrons 
in  the  epi-layer  region  outside  the  space  charge  region.  Ho- 


74 
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Fig.  5.  Correction  factor  7C  [6]  for  several  anode  radii  and  doping 
densities.  Other  parameters  for  this  diode  are  from  table  I D1038. 
The  ohmic  contact  is  at  the  rear  of  the  device  with  a radius  of 
rsc  + 4 fim. 

wever,  this  has  an  adverse  effect  on  the  breakdown  voltage. 
The  impact  of  saturation  velocity  on  the  performance  of 
diode  UVa  — 6P4  from  the  University  of  Virginia,  table  I, 
has  been  demonstrated  in  [4]  at  a fundamental  frequency 
of  100  GHz . This  diode  suffers  from  velocity  saturation  at 
this  frequency  because  of  its  low  doping  concentration  in 
the  epi-layer. 

This  effect  is  also  observed  in  highly  doped  devices  at 
higher  frequencies  as  illustrated  in  figure  6,  where  the  edge 
of  the  depletion  region  can  follow  the  voltage  waveform  at 
100  GHz,  but  cannot  follow  the  RF  voltage  waveform  at 
700  GHz.  The  diode  contact  is  at  right-hand  side  in  figure 
6.  The  diode  is  similar  to  D1038  with  Rsc  — 2.5  /xra, 
but  with  different  epi-layer  doping  This  point  can  be  bet- 
ter observed  in  figure  7 through  contours  of  the  space 
charge  region  edge  for  100  GHz  and  700  GHz.  The  velo- 
city saturation  effect  produces  a smaller  modulation  of  the 
space-charge  region  during  one  period  of  the  exciting  sig- 
nal, which  in  turn  leads  to  a smaller  capacitance  modula- 
tion, with  the  corresponding  decrease  in  the  efficiency  in 
frequency  multipliers.  Therefore,  a maximum  frequency 
for  the  pump  signal  can  be  defined  at  which  velocity  satu- 
ration effects  become  important.  It  is  also  interesting  to 
point  out  that  the  modulation  of  the  space  charge  region 
has  hysteretic  behaviour  under  high  frequency  pumping  [8]. 

Figure  8 presents  the  maximum  pump  frequency  versus 
doping  concentration  calculated  with  the  analytical  model 
presented  in  [7]  and  with  the  numerical  model  introduced 
here.  The  analytical  model  is  based  on  simple  equations  for 
depletion  layer  depth,  critical  electric  field,  and  the  avalan- 
che breakdown  voltage  (the  voltage  swing  is  selected  such 
that  it  sweeps  from  the  breakdown  to  forward  conduction 
Vrf  = (Vo  — Vbd]DC ))•  In  the  numerical  model,  the  calcu- 
lated maximum  pump  frequency  is  defined  as  the  frequency 
where  the  amplitude  of  the  calculated  current  is  25  % less 
than  the  extrapolated  amplitude  for  lower  frequencies.  Gi- 
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Fig.  6.  Variation  of  the  edge  of  the  depletion  region  with  applied  RF 
voltage  at  a)  signal  frequency  of  100  GHz,  b)  signal  frequency  of 
TOO  GHz.  The  epi-layer  doping  is  Nepi  = 4 • 1017/cm3,  the  epi- 
layer  thickness  is  tepi  = 560  nm,  and  the  substrate  length  is  80 
fjLTn  (the  Schottky  contact  is  on  the  right-hand  side). 

ven  a doping  concentration  of  Nepi  = 2 * 1017/cm3  the  ma- 
ximum pump  frequency  without  velocity  saturation  effects 
is  found  to  be  fc  < 210 GHz. 

It  is  possible  to  increase  the  maximum  pump  frequency 
by  decreasing  the  magnitude  of  the  bias  voltage.  For  a 
given  doping  concentration,  the  maximum  pump  frequency 
can  be  varied  by  changing  the  bias  voltage  and  the  RF 
voltage  swing  as  can  be  demonstrated  in  figure  9 both  for 
the  analytical  and  numerical  method.  This  implies  that 
velocity  saturation  effects  can  be  avoided  by  decreasing  the 
bias  voltage  for  a given  doping  concentration. 

IV.  Conclusions 

Two-dimensional  simulator  is  a helpful  tool  for  the  de- 
sign and  optimization  of  Schottky  diodes  both  from  the 
geometrical  point  of  view  and  from  physical  considerations. 
The  two-dimensional  simulator  allows  to  minimize  parasi- 
tics  like  series  resistance  by  a proper  design  of  the  geometry. 
Besides,  a more  accurate  computation  of  the  capacitance 
taking  into  account  edge  effects  will  result  in  a more  accu- 
rate design  of  frequency  multipliers. 
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Fig.  7.  Time  waveforms  of  the  depletion  region  edge  corresponding  to 
figure  6 a)  signal  frequency  of  100  GHz,  and  b)  signal  frequency 
of  700  GHz. 

400 
350 
300 
250 
200 
150 
100 
50 

0 0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5 

Nd  x 1017(cm~3) 

Fig.  8.  Maximum  pump  frequency  versus  epi-layer  doping  concentra- 
tion for  optimum  operation. 


Fig.  9.  Maximum  pump  frequency  versus  bias  voltage  V^a3  for  three 
different  doping  profiles.  The  RF  swing  is  chosen  as  Vrf  = 2 X 
(Vb  - Vbia3) . The  continuous  lines  are  analytical  results  from 

[7]. 

The  limiting  mechanisms  for  high  frequency  operation  of 
Schottky  diodes  have  been  discussed.  It  has  been  demons- 
trated that  saturation  velocity  effects  could  be  omitted  for 
well  designed  diodes  at  all  operating  conditions. 
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